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in ! ABSTRACT 

00 . We report the discovery of the optical counterpart of the 12.3 s transient X-ray pulsar 

O ■ GS 0834-43. We re-analysed archival ROSAT PSPC observations of GS 0834-43 ob- 

taining two new refined positions, ~14" and ^18" away from the previously published 
one, and a new spin period measurement. Based on these results we carried out optical 
and infra-red (IR) follow-up observations. Within the new error circles we found a 
relatively faint (V=20.1) early type reddened star (V~R=2.24). The optical spectrum 
shows a strong Ha emission line. The IR observations of the field confirm the pres- 
Qh| ence of an IR excess for the Ha-emitting star (K'=11.4, J-K'=1.94) which is likely 

Q . surrounded by a conspicuous circumstellar envelope. Spectroscopic and photometric 

' data indicate a BO-2 V-IIIe spectral-type star located at a distance of 3-5 kpc and 

^ I confirm the Be-star/X-ray binary nature of GS 0834-43. 
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1 INTRODUCTION 10-100) and low-luminosity (<10^^ ergs~^) pulsed persis- 

„ , , , tent emission (Negueruela 1998). 4U 0115+634 (P=3.6s), 

Be/X-ray bmary systems (BeXBs) represent the majority ^ ^332^3^ ^p^^^ 2030+375 (P=41.7 s) are all 

01 the known Hish Mass X-ray Bmanes (HMXBs) host- , r n , . 1 a .11.. 

^ / examples or the hrst two classes. Among the latter group, 

me an accretmg rotatmg magnetic neutron star (White, „ , ^ 1 ^ r> oor \ 

„ there are the well known X-ray pulsators XPer (P=835sl 

Nagase & Parmar 1995). Phenomenologically, in the X- RX J0146.9+6121 (P=1455s). Moreover, differences be- 

ray band, BeXBs can be divided into at least three sub- , , , ^ 1 1 , 1 n t n . ^1 1 1 , . r 

tween the Galaxy and the Magellanic Cloud population ot 

classes: (1) bright transients which display giant X-ray out- , , ,1.1 1 1 1 1 1 ^ i-rr 

^' ° 00 „i / i- .7 o J BeXBs have been studied and are probably due to the diner- 
bursts up to Li = 10 ergs (Type II; Stella et al. 1986) _i ^ r j-- ^ / ^ 1 mnnN n • • 

, ,,,,,, ent star formation rate (see otevens et al. 1999). B-emission 

unrelated with the orbital phase, with high spm-up rates, > , , , , , , ■ , 1 1 . 1 , , ■ , 

(Be) spectral-type stars are characterized by high rotational 

(11 transients which display periodic outbursts ol rela- 1 • / , r,noy r ii • u i 1 \ j u 

^ ' , qfi IT 1 < velocities (up to 70% ol their break-up velocity), and by 

tively high luminosity ~ 10^'' -10^^ ergs" ^ Type I) gen- . , r + - i • +• 1, 1, ■ 1,/ j 

, , , episodes ol equatorial mass ejection which might produce a 

erally occurring close to the periastron passage ol the neu- ... . r 1 ^ ^ ■ ^ ^- ■ ^ 

° f- o rotating ring or gas around the star at irregular time mter- 

tron star, and (iii) sources displaying no outbursts, but t \j_ j_- , 1 n j_ ^-m ^ 1 -il 

^ ' f- J o vals. At optical wavelengths. Be stars are difficult to classify 

comparatively moderate variations (up to a factor or ~ 1 ^ j_i ■ j_ n 1 -i 1 

^ due to the presence ot the circumstellar envelope responsible 

for the emission-lines. 

The hard X-ray transient GS 0834-43 262.0, 

* Partially based on observations carried out at ESO, La Silla, 6/1 ~-1.51) was discovered by the WATCH experiment on- 

Chile (62.H-0513) board GRANAT in 1990 at a flux level of about 1 Crab in the 

t Affiliated to ICRA. 5-15 keV energy band (Sunyaev 1990). The source was later 
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observed by GINGA (Makino 1990a, 1990b) and ROSAT as 
a part of the All Sky Survey (Hasinger et al. 1990) providing 
a much better positional accuracy (radius of 50"). Moreover 
pulsations at a period of 12.3 s were observed during the 
GINGA, ROSAT and ART-P observations (Makino 1990c; 
Aoki et al. 1992; Hasinger et al. 1990; Grebenev & Sunyaev 
1990). A target of opportunity ROSAT pointing performed 
in 1991 May allowed the source position to be determined 
with an uncertainty radius of 9" (Belloni et al. 1993). An 
optical follow-up observation of the stars within this error 
circle was performed at the ESO New Technology Telescope 
(NTT) in January 1991: no plausible optical counterpart 
was found down to a limiting magnitude of R~23.5 (Belloni 
et al. 1993). 

GS 0834-43 was also monitored by the Burst And Tran- 
sient Source Experiment (BATSE) on the Compton Gamma 
Ray Observatory (CGRO) between April 1991 and July 
1998. In particular seven outbursts were observed from April 
1991 and June 1993 with a peak and intra-outburst flux 
of about 300mCrab and <10mCrab, respectively (Wilson 
et al. 1997). The recurrence time of 105-115 days was in- 
terpreted as the orbital period of the system. However, no 
further outbursts have been observed since July 1993 either 
with CGRO/BATSE and the All Sky Monitor (ASM) on 
board the Rossi X-ray Timing Explorer. All these findings 
suggested that GS 0834-43 is a new Bc-star/X-ray binary 
system in an eccentric orbit (Wilson et al. 1997). However 
the lack of any plausible optical counterpart remained a key 
point against this interpretation. 

In this paper we report the discovery of the optical 
counterpart of GS 0834-43 a V—20.4 Be spectral-type star. 
Starting from three public ROSAT PSPC observations of 
GS 0834-43 we obtained two new independent position mear- 
surements (uncertainty radius of 10"). Within these new 
error circles we found a highly reddened star (V-R=2.24), 
the optical spectrum of which shows a strong Ha emission 
line. The optical counterpart was also observed in the IR 
showing that this star is by far the brightest object in the 
field (if'=11.4). Our findings together, obtained from ob- 
servations in three different energy bands, support the Be- 
star/X-ray binary nature of GS 0834-43. 



2 X-RAY OBSERVATIONS 

The PSPC (0.1-2.4 keV) detector on board ROSAT ob- 
served the field including GS 0834-43 several times. In the 
ROSAT public archive there are 5 observations performed 
between April and December 1991. PSPC images were ac- 
cumulated in the 0.5-2.0 keV range in order to reduce the 
strong and spatially inhomogeneous local background dom- 
inated by the soft X-ray emission from the Vela Supernova 
Remnant. 

Both a sliding cell and a Wavelet transform based de- 
tection algorithm were used in order to characterize the 
physical parameters (position, count rate, S/N ratio, etc.) 
of GS 0834-43 when detected, and to obtain a 3a count rate 
upper limit in case of non-detection (Lazzati et al. 1999; 
Campana et al. 1999). Table 1 summarizes the results of this 
analysis. 

GS 0834-43 was detected in May and December 1991 
observations as a highly variable source: from 0.067 ct s~^ 
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Figure 1. The 1991 December 17-21 ROSAT PSPC light curve 
folded at the best period P=12.307s. 

(May 5) up to 1.14 ct s~^ (December 21) corresponding 
to a variation of a factor of ~20. We note that the flux we 
obtained for sequence 160062 is consistent with that inferred 
by Belloni et al. (1993). The source was also detected during 
sequence 160061 (see Tablet). 

For each ROSAT observability window of GS 0834-43 
we obtained an independent source position. 
These were determined to be R.A.=08'' 35"" 55!6, Dec.=- 
43° 11'07'.'3 (sequence 160061; May 91; equinox 2000), and 
R.A.=08'' 35"* 55!2, Dec.= -43° 11' 10'.'3 (sequence 500127- 
8; Dec 17-21; equinox 2000). The statistical uncertainty cor- 
responds to an error radius of only 0"5 and 0'.'2 for May 5 
and December 17-21, respectively (90% confidence level). 
However due to the uncertainty in the borcsight correction 
the error radius increases to 10". We excluded from our po- 
sition analysis sequence 160062 as the source was close to 
the circular support rib (radius of ~20') of the PSPC, re- 
sulting in a higher uncertainty (error radius of ^20"). We 
note that the new refined positions of GS 0834-43 lie about 
-18" (160061) and ~15" (500127-8) away from that ob- 
tained from sequence 160062 (Belloni et al. 1993; error ra- 
dius 9"). 

The ROSAT event list and spectra of GS 0834-43 were 
extracted from a circle around the best X-ray position (with 
an extraction radius corresponding to an encircled energy of 
90% at the relevant off'-axis angle). The photon arrival times 
were corrected to the baryccntre of the solar system and 
a background corrected 1 s binned light curve accumulated 
for each observation. In order to maximize the period search 
sensitivity we analyzed only the observation with the highest 
statistics, namely sequence 500127-8 merged. A power spec- 
trum was calculated over the entire observation duration 
following the method described by Israel & Stella (1996). 
To increase the search sensitivity we searched for significant 
peaks in a narrow period interval centered around the pe- 
riod value detected by BATSE between 1991 December 28 
and 1992 Jarmary 2 (Wilson ct al. 1997; 12.307s), thus only 
a few days after the ROSAT observation. We found a sig- 
nificant peak at a confidence level of 99.5% at a period of 
12.307±0.003 s (90% uncertainties are used through out this 
paper). The pulsed fraction is 7%±2%, while the 0.1-2.4 
keV pulse shape is well fitted by two sinusoidal waves (see 
Fig. 1). A similar shape and pulsed fraction value were ob- 
tained at higher energies in 1991 December by BATSE (see 
Fig. 5, 11 and 12 in Wilson et al. 1997). 
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Table 1. ROSAT and ASCA observations of GS 0834-43. 



Pointing 

iNTiTm rifiT 


Instr. 


Expos. 


Start Time 


Stop Time 


Count ratet 
(ct s"-') 


Off-axis 


Coordinates 


500015 


PSPC 


1854 


1991 Apr 25 13:05 


1991 Apr 25 13:30 


<0.20 


44 




160062 


PSPC 


2191 


1991 May 05 16:41 


1991 May 05 17:18 


0.510±0.025 


22 


see Belloni et al. 1993 


160061 


PSPC 


2005 


1991 May 05 19:46 


1991 May 05 20:21 


0.067±0.018 





R.A.=08'' 35™ 55!6 
















Dee.=-43° 11'07'.'3 


500015 


PSPC 


8666 


1991 May 06 06:40 


1991 May 06 14:58 


<0.11 


44 




500128 


PSPC 


1316 


1991 Die 04 01:19 


1991 Die 04 01:41 


0.237±0.023 


12 




500127-8 


PSPC 


18235 


1991 Die 17 11:31 


1991 Die 21 08:38 


1.137±0.019 


12 


R.A.=08'' 35™ 55?2 
















Dee.= -43° 11' 10'.'3 


42023000 


CIS 


35916 


1994 Nov 28 07:45 


1994 Nov 28 19:01 


<0.01 


10 





t ROSAT PSPC and ASCA CIS loeally background correeted count rates were obtained in the 0.5-2.0 keV and 0.5-10 keV energy 
bands, respectively. The ROSAT PSPC count rates are vignetting and PSF corrected. Errors are at Icr level while upper limits 
at 3iT. 

t Equinox 2000. Error radius of 10" at 90% confidence level. 



Table 2. ROSAT PSPC spectral results for GS 0834-43 



Parameter 


May 1991 


May 1991 


Dec 1991 




(160062) 


(160061) 


(500127-8) 





Power-law model 




Nh (10^2 em-2) 






2-3±8:i 


r 


0.7 (fixed) 


0.7 (fixed) 


0.7±0.7 


Fx (10-11) 


1.73 


0.17 


5.46 


Lx (Nh=0; 10^5) 


2.58 


0.13 


8.05 




Blackbody model 




Nh (1022 cm-2) 




o.8±»:| 


2.o±E5:i 


kT (kcV) 


l.O(fixed) 


l.O(fixcd) 




Fx (10-11) 


1.71 


0.16 


5.44 


Lx (Nh=0; 1035) 


1.781 


0.10 


5.50 


Radius (km @ 5kpe) 


9.2 


2.1 


16.7 



Note. — The X-ray fluxes (units arc erg cm-2 s-i) ^j^g 
absorbed luminosities (units are crgs-i; a distance of 5kpe was 
assumed) refer to the 0.1-2.4 keV energy band. 



The PSPC Pulse Height Analyser (PHA) rates were 
grouped so as to contain a minimum of 20 photons per 
energy bin. The spectra of sequences 160062, 160061 and 
500127-8 (December 17-21) were fitted simultaneously. By 
using the brightest spectrum parameters as templates we 
fitted the three datasets keeping fixed the photon index and 
the temperature of the power-law and blackbody, respec- 
tively, for the two fainter spectra. A power-law (r=0.7) 
as well as a blackbody (kT=l keV) model gave good fits 
(xJ=1.09 and y^ = \.ll. respectively). 

The field of GS 0834-43 was also observed on 1994 
November 28 with the ASCA satellite. The source was not 
detected in the 0.5-10 keV band and a 3cr upper limit of 
0.01 ct s~i was inferred (corresponding to ~4xl0~i3erg s~i 
cm~^ assuming the best power-law model parameters of Tar 
ble2). 



3 OPTICAL FOLLOW-UP 

Based on the new X-ray positions, we carried out an optical 
follow-up from the ESQ (La Silla, Chile) on 1999 January 
18-20 with the 1.5 m Danish telescope and on March 13 with 
the New Technology Telescope (NTT). 

Imaging and photometry in V, R and Gunn-i bands 
(200 s exposure time each) were performed on 1999 January 
18-19 with the Danish Faint Object Spectrometer Camera 
(DFOSC), while spectroscopy of the brightest stars within 
the X ray position uncertainty regions was obtained with 
the same instrument on 1999 January 19-20. The data were 
reduced using standard ESO-MIDAS procedures for bias 
subtraction, flat-field correction, aperture photometry and 
one dimensional stellar and sky spectra extraction. Profile 
fitting photometry was also carried out with DAOPHOTII 
(Stetson 1987). Cosmic rays were removed from each frame, 
and the spectrum corrected for the atmospheric extinction 
and flux cahbrated. Figure 2 shows the fleld of GS 0834-43 
in the R filter (left panel) with the new X-ray uncertainty 
circles obtained from the 1991 May and December ROSAT 
observations superimposed. The earlier X-ray error circle 
(labelled as OLD) is also shown. Stars A and star C lie out- 
side the new uncertainty regions (~2" and -^6" away, respec- 
tively) while star D is the only object consistent with both 
circles and it is located at R.A. = 08** 35™ 55';4 and Dec. = - 
43° 11' 11'.'9 (equinox 2000; cahbrated with DSSl plates; un- 
certainty 1"). Photometry measurements for stars A, B, C 
and D are reported in Table 3. 

Since the spectroscopic properties of star A were in- 
vestigated by Belloni et al. 1993 (star A is a late spec- 
tral type star without amy significant emission features), we 
focused our attention on star D. We obtained three low- 
resolution (llA) spectra (30min exposure time each; l'.'5 
slit; l'.'2 seeing) of star D on 1999 January 20 with a grism 
covering the spectral range 5200-lOOOOA (see upper line of 
Fig. 3) . Note that the blue part of the spectrum, being rel- 
atively faint (V=20.4), was below the instrument sensitivity 
(telescope-|-CCD-|-grism). After reduction, the spectra were 
summed to increase the S/N ratio. Star D was also observed 
on 1999 March 13 with the ESO Multi-Mode Instrument 
(EMMI) mounted on the adaptor-rotator at the Nasmyth 
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Figure 2. R-band (left) and K'-band (right) images of the field of GS 0834-43 together with the new X-ray position uncertainty 
regions (10" radius) obtained with the 1991 ROSAT PSPC observations. North is top, east is left. Stars D represents the proposed 
optical counterpart of GS 0834—43. 



Table 3. Optical and IR results for star A,B,C and D 



Star 


V 


R 


V-R 


J 


H 


K' 


J-K' 


A 


18.9 


17.8 


1.1 


15.86 


15.06 


14.80 


1.06 


B 


20.4 


19.0 


1.4 


16.20 


15.24 


14.98 


1.22 


C 


22.9 


20.8 


2.1 


15.40 


13.81 


13.29 


2.11 


D 


20.4 


18.2 


2.2 


13.33 


12.26 


11.39 


1.94 



Note — Optical magnitude absolute uncertainty is 
~0.02. 



-0.2, IR is 



B focus of the NTT. A 30min low-resolution (llA; 2" slit; 
2" seeing) spectrum was obtained in the 4000-8500A range 
(see Fig. 3; lower line). 

Due to its faintness, star D is difficult to study. The 
S/N ratio of the summed DFOSC spectrum is acceptable 
only above ~6000A, while at >7500A the signal is domi- 
nated by interference fringes. The steep rise of the spectrum 
in the UV argues for a very hot and reddened object. The 
absence of strong forbidden lines rules out the possibility of 
a pre-main sequence object or a cataclysmic variable. More- 
over the absence of strong absorption features points to an 
OB spectral-type star. The main features of the DFOSC 
spectrum are: (i) a strong Ha emission-line with equiva- 
lent width EW~-30±2A, and (ii) pronounced interstellar 
absorption hues for Nail (5890A and 6270A). The Ha line 
shows an extended and asymmetric profile. At llA spectral 
resolution, this might indicate a line splitting, possibly due 
to the presence of a disk. All this evidence clearly points to 
the association of this object with the X-ray source. 

The EMMI spectrum has a slightly better S/N ratio 
and allows us to extend the spectral region into the 5500- 
8400A range. However, a precise spectral classification is 
not yet possible. A few Sill lines (e.g. the strong 5466.43- 



5466. 87A doublet) are recognizable, suggesting a B2 spectral 
type. However the uncertainty on the spectral classification 
is large, ranging from an 08e to a B3e star. Emission features 
corresponding to He I lines (e.g. 6678 and 7281 A) are also 
evident, while the He I lines at 7065 and 7818A look filled in, 
as well as faint nebular [N H] lines on both sides of Ha, which 
most likely are due to the superposed emission from the 
Vela SNR (see Fig. 3). In the EMMI spectrum the Ha EW 
is -33±2A, while the FWZH is 41 ±3 A, corresponding to 
a stellar wind terminal velocity Y^ind —1800 km s~^. Deep 
interstellar lines and bands are clearly visible in the star 
spectra, the NaD2 doublet having an EW of 4.1±0.8A. 

Star B the position of which is still marginally consis- 
tent with the new X-ray error circle was observed on 1999 
January 19. We obtained a low-resolution (llA) spectrum 
(1 h exposure time; 2'.'0 slit; l'.'7 seeing) with a grism covering 
the 3500-7000A spectral range. Star B is even more difficult 
to study being fainter than star D. Although the spectrum 
has a low S/N ratio, we detected large and deep absorption 
features typical of late type stars (probably an early M). 
Moreover no emission lines were visible in its spectrum. 



4 IR OBSERVATIONS 

Infrared observations were carried out on 1993 April 1 dur- 
ing service observing time using the Infrared Imager Spec- 
trograph (IRIS) instrument (Allen 1993) mounted at the 
Cassegrain focus of the 3.9m Anglo-Australian Telescope 
(AAT), Siding Springs Observatory, New South Wales. Im- 
ages were obtained in the J (60s), H (20s) and K' (5 sees.) 
bands. The data were reduced using the Starlink Figaro 
software (Shortridge et al. 1998), and analysed using the 
PHOTOM package (Eaton & Draper 1998). The standard 
stars HD 84090 and HD 100231 were used for calibration 
purposes, giving absolute uncertainties of around 0.02 mag- 
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resolution spectra of GS 0834 43 obtained on 1999 January 20 (upper line; Danish+DFOSC) and March 13 (lower line; 
Fluxes are arbitrarly normalised and a constant was added to the January 20 one. 



nitudes in each band (although the H band image was near 
the saturation limit of the detector and should be treated 
with a little more suspicion). 

Taking the colours of a B2V star from Wegner 
(1994) as (V-K)=-0.67, (J-K)=-0.18 and (H-K)=-0.04, us- 
ing the reddening relationships which gives E(J-K)/E(B- 
V)=+0.50, and combining these with the observed values for 
GS 0834-43, we find E(J-K)=+2.0 and thus E(B-V)=+4.0. 
Taking the parameter Rea!t=3.3 (Rea!t=Av/E(B-V)), we es- 
timate an extinction of Av=13.2 from the IR magnitudes. 



5 DISCUSSION 

The X-ray, optical and IR observations of the field of 
GS 0834-43's presented here, led to the identification of 
the optical counterpart of this 12 s transient X-ray pulsar 
discovered in 1990. The measurement of the distance to 
GS 0834-43 based on the optical data is hampered by the 
uncertainties in the spectral classification. 

However some information can be inferred based on 
our optical and IR photometric measurements. The intrin- 
sic 1^ — i? color for GS 0834-43 is ~ —0.1 (assuming a main 
sequence star with spectral type in the 09-B3 range). Since 
the observed V—Ris^ 2.25 the reddening should amount to 
_Ev-R ~ 2.35, and assuming a standard reddening law (Fitz- 
patrick 1999) this converts to Ar ^ 7 (regardless of whether 
the reddening medium is uniformly distributed along the line 
of sight or is intrinsic to the source). Moreover from the X- 
ray spectral fits, we derive a Nh of 1-2 x 10^^ which is 
in good agreement with the photometric reddening estimate 
(Bohhn et al. 1978). The observed magnitude of GS 0834-43 
is J? ~ 18.2 and the absolute R magnitude of a star with the 
assumed spectral type is within —2.5and — 5. Taking into 
account the effects of reddening, a distance modulus within 
13.7 — 16.2 can be derived. This translates to a distance of 
~5-17kpc, well beyond the Vela SNR located at ~500pc. 
However the Galaxy edge in the direction of GS 0834-43 is 
located at ~6kpc, so we can put a limit on the spectral- 



type of star D assuming that it is placed at the Galactic 
border. For an 09, BO or B2 dwarf star we obtain a red- 
dened {Ar ^ 7) R magnitude at 6 kpc (distance modulus 
of ~ 14) of 16.5, 17 and 18.5, respectively. When compared 
with the observed R magnitude of 18.2 for star D, this sug- 
gests a B2Ve spectral-type star. 

Also the derived E(B-V)=-l-4.0 from the IR data sug- 
gests a highly reddened object. Using the relationship for the 
contribution of the circumstellar (cs) environment derived 
by Fabregat and Reglero (1990), where E^s^O. 0049-0.00185 
EW(Ha), and given the observed EW(Ha)'— 30A, we es- 
timate Ecs to be only 0.06 magnitudes. This suggests that 
the bulk of the extinction is interstellar rather than local 
to GS 0834-43. This results confirm that the contribution of 
the disc (which is clearly present, given the large equivalent 
width inferred for Ha) to the observed IR colors is relatively 
small. 

However the value of Av^lS derived from the IR ob- 
servations differs significantly from that estimated from the 
optical data (where A_fj~7, and thus A\/~9.4), suggesting 
that the object is either earlier in spectral class than B2, 
closer that 5-6 kpc, or possibly not a main sequence ob- 
ject. The observed magnitudes are similar to those of EXO 
20304-375 (see e.g. Negueruela 1998), which is compatible 
with a BOV star at ~3 kpc or a BOIII at ~5 kpc. 

Based on both photometric and spectroscopic findings 
we conclude that star D is most likely a BO-2 V-IIIe star at 
a distance of 3-5 kpc. A more accurate distance and spec- 
tral clEissification will have to await detailed optical spectro- 
scopic observations in the blue end of the spectrum. 

For a distance of 5 kpc and extrapolating the 0.1- 
2keV ROSAT fiuxes to the 1-10 keV band (we as- 
sume the spectral model fitted by Aoki et al. 1992 using 
Ginga data) we obtain an outburst peak luminosity Lx (1- 
lOkeV) ~ 5 8 X 10'"' erg s^"'. A somewhat higher luminosity 
(~10'^^erg s^^) would be inferred by extrapolating the 20- 
100 keV BATSE peak flux reported by Wilson et al. (1997) 
to the 1-10 keV band. Such a peak luminosity is a typical 
value shown by X-ray pulsars which also display Type I out- 
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bursts (Stella et al. 1986; Negueruela 1998) occurring close 
to the time of periastron passage and with a periodic recur- 
rence given by the orbital period of the system. 
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